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The spectral weight functions and the optical conductivity of the Holstein model are studied on 
a one-dimensional six-site lattice with periodic boundary conditions for three different electron 
concentrations: a single electron, two electrons of opposite spins, and half filling. A density matrix 
approach is used to obtain an optimal phonon basis and to truncate the phonon Hilbert space 
without significant loss of accuracy. This approach allows us to calculate spectral functions for 
electrons dressed locally by the optimal phonons as well as for bare electrons. We obtain evidence 
for a smooth crossover from quasi-free electrons to an heavy itinerant small polaron (single-electron 
case) or bipolaron (two-electron case) as the electron-phonon coupling strength increases. At half 
filling we observe a crossover from a quasi-free-electron ground state to a quasi-degenerate Peierls 
charge-density-wave ground state for a finite electron-phonon coupling. This crossover is marked by 
an abrupt drop of the Drude weight which is vanishingly small in the Peierls phase. 

PACS Numbers: 71.38. +i, 71.10.Fd, 71.10.Pm, 63.20.Kr 



I. INTRODUCTION 

The Holstein modelS has been used for many years to 
study physical problems related to the electron-phonon 
interaction, such as the formation of polarons and bipo- 
larons by self-trapping of charge carriers, or the exis- 
tence of charge-density-wave (CDW) ground state due 
to the Peierls instability.El While our knowledge of the 
ground state of this model has considerably progressed 
for the past few years, our understanding of its dynam- 
ical properties is still very limited and often disputed. 
The lack of reliable results is especially important in 
the non-adiabatic and intermediate electron-phonon cou- 
pling regimes, where most of the interesting physics oc- 
curs, such as self-trapping crossover and metal-insulator 
transition. Currently, there is no well-controlled ana- 
lytical method to study these regimes and most reliable 
results come from-aurnerical simulations, such as exact 
diagona^ationsQLj, quantum Monte Carlo (QMC) sim- 
ulational3il3, and recent deusitaz-matrix renormalization 
group (DMRG) calculationsOEl Among these various 
numerical methods, only exact diagonalizations can eas- 
ily be used to compute dynamical properties of the sys- 
tem. Although this technique can only be applied to 
small clusters due to restrictions on computer resources, 
it often allows us to gain a useful insight into the physics 
of the system. Moreover, if the error due to the necessary 
truncation of the phonon Hilbert space is negligible, this 
method provides numerically exact results, which can be 
used to assess the accuracy of other analytical or numer- 
ical methods. 

In this paper we report our study of the dynamical 
properties of a six-site one-dimensional Holstein lattice 
with periodic boundary conditions. We consider three 
different electron concentrations: a single electron, two 
electrons of opposite spins and a half-filled band (six 
electrons with zero total spin). The Holstein model de- 
scribes non-interacting electrons coupled to dispersion- 
less phonons. Its Hamiltonian is 



h = n b h - 7 fa + 6 *) ni 

i i 

where c\ a (ci a ) creates (annihilates) an electron with spin 
a on site i, n.i — cLcjf + c^c^, and b\ and hi are 
creation and annihilation operators of the local phonon 
mode. The model parameters are the hopping integral 
t, the electron-phonon coupling constant 7, and the bare 
phonon frequency f2. For all the results presented in this 
paper, the phonon frequency is chosen to be in the non- 
adiabatic regime fl = t and we study the variations of the 
system properties when the electron-phonon coupling 7 
goes from zero to the strong-coupling regime 7 > f2, t. 

We perform exact diagonalizations of the Holstein 
Hamiltonian using the efficient local phonon Hilbert 
space |teduction method that we have recently intro- 
duced]^ This approach uses the information contained in 
a reduced density matrix to generate an optimal phonon 
basis which allows us to truncate the phonon Hilbert 
space without significant loss of accuracy. In our pre- 
vious work, this method has been demonstrated on the 
ground state of the six-site Holstein model at half filling. 
Here, we extend this approach to different band fillings 
and to excited state calculations. We also show how to 
use the optimal phonon basis to dress electrons with local 
phonons. Using the Lanczos algorithmca, single-particle 
and pair spectral functions are calculated for bare elec- 
trons and dressed electrons, and the Drude weight and 
optical conductivity are computed. 

This paper is organized as follows: in the next section, 
we describe our method to obtain an optimized phonon 
basis and to dress electrons with these phonons, and in- 
troduce the dynamical quantities we have calculated. In 
Sec. Ill we present our results for the six-site Holstein 
model. Finally, Sec. IV contains our conclusions. 
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II. METHODS 



A. Optimal phonon basis 



In order to perform an exact diagonalization of the 
Hamiltonian ([!]), one needs to introduce a finite basis 
to describe the phonon degrees of freedom. If one uses a 
bare phonon basis (the basis made from the lowest eigen- 
states of the operators &J&i), the number of phonon levels 
needed for an accurate treatment can be quite large in 
the strong-coupling regime. However, this number can 
be strongly reduced by using an optimal basis (a basis 
that minimizes the error due to the truncation of the 
phonon Hilbert space). In a previous worka we have 
introduced a density matrix approach for generating an 
optimal phonon basis. The key idea, of this approach is 
identical to the key idea of DMRQH3: in order to elimi- 
nate states from a part of a system without loss of accu- 
racy, one should transform to the basis of eigenvectors of 
the reduced density matrix, and discard states with low 
probability. To be specific, consider any wave function 
\t/j) in the Hilbert space of the Holstein model. Let a la- 
bel the four possible electronic states of a particular site 
(empty, occupied by a single electron of spin up or down, 
or occupied by two electrons of opposite spins) and let n 
label the bare phonon levels of this site. Let j label the 
combined states of all of the rest of the sites. Then \tp) 
can be written as 



(2) 



The reduced density matrix p of the state for this site 
is 



E 



|a)(a| (g) 



\n)(r\ 
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where r is another index labeling the bare phonon levels. 
This density matrix is always diagonal for the electronic 
states because of the conservation of the number of elec- 
trons. The phonon density matrix for each electronic 
state a of the site is given by 



Pn,r — ^ ; jV'ar,/ 
3 



(4) 



Let w a k be the eigenvalues and <j) a k (n) the eigenvectors of 
Pn rm where k labels the different cigcnstates for a given 
electronic state of the site. The states 



14, 



(j) ak {n)\n) , k = 1,2, 



(5) 



form a new basis of the phonon Hilbert space for each 
electronic state a of the site. The w a k are the probabili- 
ties of the site being in the state 



, k) = \a)\(j> ak ) 



(6) 



if the system is in the state \ip). These states \a, k) form 
a new basis of the Hilbert space associated with each 
site. If w a k is negligible, then the corresponding state 
\a, k) can be discarded from the basis for the site, with- 
out affecting the state If one wishes to keep a limited 
number of phonon states m for a site, then the best states 
to keep corresponds to the m eigenstates of p with largest 
eigenvalues for each electronic state of the site. The cor- 
responding phonon states \<fi a k) form an optimal phonon 
basis. In the Holstein model we have found that keep- 
ing m = 3 — 5 optimal states \4> a k) per phonon mode for 
each electronic state of the site gives results as accurate 
as with hundred or more bare phonon states per site for 
a wide range of parameters. 

Unfortunately, in order to obtain the optimal phonon 
states, we need the target state (||), for instance, the 
ground state. Usually do not know this state - we want 
the optimal states to help $et it- This problem can be cir- 
cumvented in several waysu. Here, we describe the algo- 
rithm that we have used to obtain optimal phonon bases 
for the ground state and low-lying states of the Holstein 
model. First, we calculate a large optimal phonon basis 
in a two-site Holstein system with appropriate parame- 
ters. In such a small system we can carry out calculations 
with enough bare phonon levels to render completely neg- 
ligible errors due to the truncation of the phonon Hilbert 
space. Thus, target states can be obtained directly by 
diagonalization in the bare phonon basis. Then, the op- 
timal phonon states of the two-site system are used as 
the initial basis states for calculations on larger lattices. 
The simplest of algorithms described in Ref. 13 is used 
for the six-site system. A single site (called the big site) 
contains a large number of the optimal phonon states 
obtained in the two-site systems (up to a few hundreds). 
Each other site of the lattice is allowed to have a much 
smaller number of optimal phonon states, m ~ 3 — 5, for 
each electronic state of the site. Initially these states are 
also optimal phonon states of the two-site system. The 
ground state of the Hamiltonian (|l|) is calculated in this 
reduced Hilbert space by exact diagonalization. Then, 
the density matrix (||) of the big site is diagonalized. The 
most probable m eigenstates for each electronic state of 
the big site form new optimal states which are used on all 
of the other sites for the next diagonalization. These new 
phonon states are now optimized for the six-site system 
and thus are different from the optimal states of the two- 
site system. After the first diagonalization, the new op- 
timal states of the six-site system are not very accurate. 
Thus, diagonalizations of the Hamiltonian ([!]) and of the 
density matrix are repeated until the optimal states have 
converged. In each diagonalization, the big site always 
has a large number of phonon states, so that it can gener- 
ate improved optimal states for the next iteration. After 
full convergence of the optimal phonon basis, the error 
made by using 3-5 optimal states instead of hundreds of 
bare levels is negligible: typically, the error in the ground 
state energy is smaller than 10~ 5 with 3 or more optimal 
states when only the ground state is targeted. 
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To study dynamical properties we can extend the 
above approach by targeting the ground state \ipo) as well 
as some low-lying excited states \tps),s > 0. In this case, 
the density matrix of the big site is formed by adding 
the density matrices p s of each state \ip s ), with weight- 
ing factors a s (normalized by ^2a s = 1), 



P : 



E 



a s p s 



(7) 



Thus, information from several states can be included to 
select the optimal phonon basis. The weighting factors 
allow us to vary the influence of each state \ip s ) in the 
formation of the optimal phonon basis. Not surprisingly, 
we have found that more optimal states must be kept on 
each site to reach a given accuracy when several states 
are targeted than when only the ground state is targeted. 
Obviously, these additional optimal states are necessary 
to describe accurately the excited states \ips)- However, 
they do not seem to be necessary to obtain a qualitative 
desc ripti on of dynamical properties (see the discussion in 
Sec. [ID). Therefore, we usually target only the ground 
state in our calculations. 



B. Dressing of electronic operators 

A very interesting feature of our method is that it 
provides a very simple way to dress electrons with 
phonons. Let assume that the density matrix eigenstates 
4>ak(n), k = 1,2, ... are ranked by decreasing weight w a k 
for each electronic state a. For a given index k, the 
weights w a k and eigenstates <fi a k(n) of the different elec- 
tronic states a often seem completely unrelated. How- 
ever, we can consider the relative weight 



E 



(8) 



of an optimal state \<j) a k) compared to the weight of all 
states in the corresponding optimal basis. As noted pre- 
viously w a k (and thus r a k) decreases rapidly with k. We 
have found that, for a given index k the variations of r a k 
as a function of the electronic state a are much smaller 
(at least by one order of magnitude) than the variations 
of r a k between successive values of k. Therefore, there is 
an unambiguous relation between optimal phonon states 
\4>ak) for different electronic occupations of a site a given 
by their relative weights r a k- This one-to-one mapping 
between optimal states can be used to dress electronic 
operators. 

All electronic operators can be written as the sum and 
product of operators acting on a single site, for instance 
c\ a and Ci a . Such a local operator is diagonal in the bare 
phonon basis and can be written as 




a)(P\ ®I P h 



(9) 



where a and label the four possible electronic states 
of the site and I p h is the identity operator acting in the 
Hilbert space of the local phonon mode. Now we can 
define the corresponding dressed operator as 



O 



Y,(° a ,P \a){f3\®Ua, ) , 

a,/3 



(10) 



where U a ,p is a unitary operator in the phonon Hilbert 
space given by 



U a ,l3 = ^2 l^afe)(^/3fc 



(11) 



Obviously, for a — (3 we have U a ,p — Iph because the 
density matrix eigenstates 4> a k(n) satisfy the orthonor- 
malization condition 4>ak(n) 4> a k{r) = b~n,r- However, 
for a ^ (3, U at p is not trivial because the eigenstates for 
different electronic states are unrelated, in general. Nev- 
ertheless, U a ,f3 is unambiguously defined thanks to the 
one-to-one mapping between eigenstates discussed above. 
To be rigorous, the operator U a .p is unitary only if we 
sum up the index k in Eq. ( [Tl| ) over the infinite num- 
ber of states in the basis (|^) . As we generally know only 
to ~ 3 — 5 optimal phonon states, the operator U a ,p also 
involves a projection onto the subspace spanned by these 
few states. However, by construction all wave functions 
that we calculate have a negligible weight out of this sub- 
space. Thus, U a ,p can be regarded as a unitary transfor- 
mation for all practical purposes. 

Clearly, the operator O transforms electronic states as 
the bare operator does, but it also transforms the phonon 
degrees of freedom accordingly. For instance, the oper- 
ator Ci iCr not only removes an electron with spin a form 
the site i, but it also transforms the phonon mode on 
this site, changing optimal phonon states for a site with 
two electrons into optimal states for a site with an elec- 
tron of spin —a and changing optimal states for a site 
with one electron of spin a into the optimal states of 
an unoccupied site. Therefore, the operator O acts on 
electrons dressed by the local optimal phonon states as 
the bare operator acts on bare electrons. For instance, 
the operator c] a creates a dressed electron of spin a on 
the site i. We note, however, that the dressing of elec- 
trons by phonons at a finite distance from the electrons 
is completely neglected with this method. 

In two cases (7 = and t = 0) it is possible to cal- 
culate the optimal phonon basis analytically and thus to 
understand the transformation (0). In the weak cou- 
pling limit (7/^,7/0 — > 0) the optimal phonon states 
resemble the bare phonon levels for each electronic state 
of the site. Therefore, the unitary transformation U a ,f3 
is similar to the identity operator for any values of a 
and (3 and thus, O ~ O. In the strong-coupling anti- 
adiabatic limit (7, >> t) the optimal phonon states 
are simply the eigenstates of quantum oscillators with an 
equilibrium position shifted by 2 r y/QN a , where N a is the 
number of electrons on the site in the electronic state a, 
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Hj I a) = N a \a). This corresponds to the states obtained 
by applying the Lang-Firsov unitary transformation^ 



S(g) = e - 9 ^ {b '- b * )n ' 



(12) 



with g = to the bare states. Therefore, in this limit 
we have 



6 = s(-f/n) o s-\j/n) 



(13) 



and the electronic operators obtained with the transfor- 
mation (|l^)|-M«|*«mpletely equivalent to those defined in 
other worksuilj'Ea using the Lang-Firsov transformation. 
However, in the general case, the transformation ( |l0|) is 
more accurate than the Lang-Firsov transformation. The 
later is only an (analytical) approximation of the trans- 
formation of a local phonon mode as a function of the 
electronic occupation of a site, while the former is based 
on a (numerically) exact transformation of the phonon 
states. 



C. Dynamical quantities 

We compute dynamical properties such as spectral 
weight functions and optical conductivity using the Lanc- 
zos algorithm combined with the continued fraction 
method.O This algorithm yields not only the dynami- 
cal correlation and response functions of the system, but 
also the most important eigenstates \tp s ) that contribute 
to these functions and that we need to build the density 
matrix (Q) . We define the spectral weight function as 



A(p, uj) = —Im 

7T 



H 



En 



-pa 



(14) 



where \ipo) is the ground state wave function for a given 
number of electrons, Eq is the ground state energy, and 
the operators c po . and c pa create and annihilate an elec- 
tron with momentum p and spin er, respectively. In all 
the results presented in this paper, we have used a broad- 
ening factor e = 0.1. The total weight of this spectral 
function (obtained by integrating over uj) is equal to the 
momentum density distribution 



n<r(p) = (^olcLcpirllpo 



(15) 



We also define a spectral function A(p,uj) and its total 
weight n a (p), where we substitute dressed operators c p(T 
and c pcr for the corresponding bare operators in ( |l4| ) and 
(|l5|). When there are more than one electron in the lat- 
tice, one expects electrons to form tightly bound pairs at 
strong electron-phonon coupling. Thus, it is interesting 
to study the pair spectral function 



P(p, uj) — —Im 

7T 
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where 



AI, 



(17) 



and its hermitian conjugate A p are pair operators with 
momentum p, and N being the number of sites. The 
total weight of this spectral function is given by 



d{p) = (?/>o|a;a p iv>o>. 



(18) 



In this case too, we introduce a spectral function P(p, uj) 
and its total weight d(p) using dressed electronic opera- 
tors At and A p in Eqs. (|l|) and (|l|). 

The real part of the optical conductivity is made up 
of a Drude peak at to = and an incoherent part for 
uj > 0, <j{uj) = D5(ui) + (t'(lu). The incoherent part of 
the conductivity is given by 



ujN 



-Im 



i 



H — Eq — uj — ie 
where the current operator is defined as 



(19) 



(20) 



<j'{uj) can be calculated using the Lanczos mcthod,E3 then 
the Drude weight D can be evaluated using the well- 
known sum rule 



f°° ire 2 
J o a(uj)duj = — (-T) 



(21) 



relating the total weight of the optical conductivity to 
the electronic kinetic energy per site 



(22) 



In this paper, the spectral functions and the incoherent 
part of the conductivity are always shown in arbitrary 
units. Quantitative results for the Drude weight D are 
expressed in units of 2ire 2 t. 



D. Comparison with other approaches 

As a first check of our method, we have compared our 
exact diagonalization results for the lowest eigenstates 
of small Holstcin clusters (up to six sites) with DMRG 
calculations. We have always found a good quantitative 
agreement. For instance, the eigenenergies obtained with 
both methods agree within 10 _3 i for at least up to the 
18th lowest eigenstates in a six-site lattice. 

We have also carried out calculations of dynamical 
properties in the two-site Holstein model where we can 
keep enough phonon levels to obtain numerically exact re- 
sults. In this case, we have calculated the optimal phonon 
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basis in various ways, changing both the number of tar- 
geted states \ip s ) and the weighting factors a s to build the 
density matrix (Q). We have found that overlaps between 
corresponding optimal phonon states in the different ba- 
sis are always larger than 90%. The spectral functions 
and conductivity calculated using the different basis also 
agree qualitatively. Therefore, the inclusion of excited 
states in the density matrix (]?]) does not seem to be nec- 
essary to obtain a qualitative description of dynamical 
properties. 

Finally we have also calculated the optical conductivity 
of the six-site Holstein model of spinless fermions at half 
filling. We have found a satisfactory agreement w-ith the 
results obtained recently by Weisse and FehskeQ using 
completely different approaches to truncate the phonon 
Hilbert space and to calculate the conductivity. 

The optimal phonon approach used in this work is so 
efficient that we can easily carry out calculations for dy- 
namical quantities that require powerful parallel com- 
puters when a standard phonon Hilbert space truncation 
method is usedjj LtEj All results presented here have been 
obtained on a workstation with a 133MHz processor and 
150Mb of RAM memory. 

III. RESULTS 
A. Single electron 

The case of the Holstein model with a single electron is 
known as the polaron problem. This case has been exten- 
sively n stjud ied with ho fk, analytical and numerical meth- 
ods. raraoHBEiyy For weak coupling (7/O < 1 
and 7 2 /J7 < 2t) the ground state is a quasi- free electron 
dragging a phonon cloud. For strong coupling > 1 

and 7 2 /f2 > 2t) the electron becomes trapped by the lat- 
tice distortion that it generates. The quasi-particle com- 
posed of this self-trapped electron and the accompanying 
lattice distortion is called a polaron. The polaron is said 
to be small when the spatial extension of the self-trapped 
state is limited to one site. It is known that a smooth 
crossover occurs from the quasi-free electron ground state 
to a small polaron ground state as the electron- phonon 
coupling increases. 

First, we examine the optimal phonon states obtained 
with our method. Figure ^ shows the optimal phonon 
wave functions 4>{q) as a function of the phonon coordi- 
nate q = b + b^ for different electron-phonon couplings. 
Only the most important optimal state is shown for 
each of the two possible electronic occupations of a site 
(N a = 0,1). For weak coupling (7 = 0.3t) the optimal 
states are similar to the bare phonon levels and thus, the 
wave function <f>(q) is just the ground state of a quantum 
oscillator. As the coupling increases, the optimal states 
change smoothly and become increasingly distinct. For 
all coupling each wave function (j)(q) has a large overlap 
with the lowest eigenstate of a quantum oscillator with a 
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FIG. 1. Single-electron system: optimal phonon wave func- 
tions as a function of the oscillator coordinate q and 
their weights W for the two possible occupations of a site 
(N a — 0, 1) and different electron-phonon couplings 7. 

shifted equilibrium position. This shift is always very 
small for N a = 0, but for N a — 1 it increases with in- 
creasing coupling and tends to 27/f} at strong coupling. 
This is in agreement with the strong-coupling theory that 
predicts optimal phonon states given by the Lang-Firsov 
transformation ( |l2| ) with g = 7/fi. However, in the gen- 
eral case, the optimal phonon states are different from 
the states obtained with this transformation. First, one 
can see in Fig. pi that the oscillator shift for N a = 1 is 
smaller than 27/ f2 for intermediate couplings (7 = Lit) 
and reaches this value only for strong coupling (7 = 3t). 
Of course, this difference could be taken into account sim- 
ply by using an effective parameter g smaller than 7/O in 
Eq. (1T3) . However, there are other features of the opti- 
mal states that a simple Lang-Firsov transformation can 
not reproduce. For instance, for 7 = 2t the optimal wave 
function 4>(q) for N a = 1 has an important tail at low q. 
This can be understood as a retardation effect due to the 
finite phonon frequency £l/t. Most of the time a site is 
unoccupied and the phonon mode is in the optimal state 
for N a = 0. When the electron hops on this site, the 
phonon mode can not adapt instantaneously. Thus, its 
state for N a = 1 becomes a combination of the states ob- 
tained using the Lang-Firsov transformation for N a = 
and N a = 1. In Fig. |l| we also give the weight Wq and 
W\ of the most important optimal state for N a = and 
N a = 1, respectively. W\ is much smaller than Wq be- 
cause the probability of finding the electron on a given 
site is only 1/6 while the probability of a site being empty 
is 5/6. For 7 = and 7 — * 00 one can show that there is 
only one optimal state for each electronic state of a site 
and thus, Wq = 5/6 and W\ = 1/6. For intermediate 
couplings, Wo and W\ become smaller, showing the in- 
creasing importance of the higher optimal phonon states. 
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FIG. 2. Single-electron system: (a) spectral functions 
A(p = 0,ui) (thick line) and A(p — 0,ui) (thin line) for different 
electron-phonon couplings 7. For 7 = 2t P(p = 0,ui) has been 
expanded by a factor 5. (b) Total spectral weight rv(p = 0) 
(open circle) and n a (p = 0) (filled circle) as a function of the 
electron-phonon coupling. 

However, for all cases presented in Fig. [I], the two most 
important optimal states constitute more than 98% of 
the total weight. 

For weak electron-phonon coupling both spectral func- 
tions A(p = 0, to) and A(p = 0, ui) have a single peak at 
lu = It because the ground state wave function is simply 

l^o) « 4=o,«r|0> « 4=0,^10), (23) 

where |0) is the vacuum state (without electron nor 
phonon), with a ground state energy of about —It. Fig- 
ure ||(a) shows that satellite peaks appear above the dom- 
inant peak energy in both spectral functions for larger 
couplings. The position of the dominant peak shift to 
higher energy as the coupling 7 increases. The distance 
between these peaks is roughly J7, with some peaks too 
small to be seen. We can easily understand the structure 
of these spectral functions. Equation ( p"4| ) with e — > 
can be written as 



A( p ,uj) =^|(0„|c pCT ^o>| 2 ^- (e n -Eo)), (24) 

n 

where Eq and |^>o) are the ground state energy and wave 
function for a single electron, and e n and \4>n) are a com- 
plete set of energies and eigenstates for a lattice con- 
taining no electron and thus, only non-interacting local 
phonons. In this case, all eigenenergies are of the type 
e n = mQ, where m is an positive integer number. There- 
fore, A(p,u>) and A(p, oS) contain only peaks with spac- 
ing O starting from — Eq. In A(p = 0,u>) the weight 
of the dominant peak shifts increasingly to the satellite 
peaks and in the strong coupling regime (7 = 2t), no 
dominant peak can be identified [Fig. 0(a)] . Moreover, 
the total weight n a (p = 0) of this spectral function de- 
creases continuously as 7 increases [see Fig. ||(b)] and in 
the strong-coupling limit it tends to the value 1/6 sig- 
naling a completely localized electron. (More precisely, 
if the electron is localized on a single site of the six-site 
lattice then n a (p) = 1/N —1/6 for all momentum p.) 
This decrease and the considerable incoherent contribu- 
tion to A(p — 0, oS) shows that the free electron state 
c p=a cr|0) i s n0 longer a good starting point for a descrip- 
tion of the ground state for 7 > 1.5t. On the other hand, 
A(p = 0, to) contains a well defined quasi-particle peak for 
all couplings [Fig. f|(a)], indicating that the motion of the 
electron is closely accompanied by a local phonon cloud 
represented by the optimal phonon states. The position 
of this peak is determined by the ground state energy Eg 
as explained above. For 7 = 2t this position gives a po- 
laron energy that approaches the strong-coupling result 
Eq = —j 2 /il. The total weight n a (p — 0) of this spectral 
function is very interesting [Fig. 0(b)]. First, n a (j> = 0) 
is always larger than 0.9, showing that the ground state 
is very well described by 

for all couplings. Furthermore, we note that n a (j> = 0) 
first decreases slightly as 7 increases and then tends to 
1 in the strong-coupling limit. We think that the ini- 
tial decrease of n a (p — 0) shows the increasing impor- 
tance of the extended phonon cloud following the electron 
as one goes from the non-interacting limit to the inter- 
mediate electron-phonon coupling regime. The operator 
Cp = <r dressed only by local phonons does not describe 
this extended phonon cloud. For larger coupling 7 the 
phonon cloud collapses to a single site as a small polaron 
is formed and thus, the dressed operator Cp =0 a becomes 
again an almost exact description of the ground state. 
In this regime, the operators & p a and & p CT obtained with 
the transformation (J10|) are similar to the "small polaron 
ope^atopj' obtained using the Lang-Firsov transforma- 
tionlf EjO as discussed in the previous section. 

The evolution of the Drude D weight and of the kinetic 
energy per site T is shown as a function of 7 in Fig. [|(a). 
The kinetic energy gives the total weight of the optical 
conductivity according to Eq. (pll) while the Drude wei- 
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FIG. 3. Single-electron system: (a) the Drude weight D 
and the kinetic energy per site T as a function of the elec- 
tron-phonon coupling 7. (b) The incoherent part of the 
optical conductivity cr'(u>) in the quasi-free electron regime 
(7 = 0.8£) and in the small polaron regime (7 = 2.5t). 

ght measures the contribution of coherent motion to the 
optical conductivity. In Fig. ||(a) the units have been cho- 
sen so that both quantities appear equal when there is no 
contribution from the incoherent part of the conductiv- 
ity er'(ii-), as in a non- interacting system (7 = 0). We see 
that both the Drude weight and the average kinetic en- 
ergy decreases smoothly as the coupling increases. How- 
ever, the Drude weight decreases much faster and be- 
comes very small for 7 > 2t. The slow decrease of D for 
small 7 shows the slight increase of the electron effective 
mass as it has to drag an increasingly important phonon 
cloud. The small but finite value of D for large 7 re- 
flects the fact that a polaron moves coherently, as a free 
carrier, in the Holstein model, but its effective mass is 
much larger than the bare mass of the electron. The de- 
crease of the ratio D/T implies that incoherent processes 
become more important in the optical conductivity as 7 
increases. As seen in Fig. |](b), er'(w) becomes also fairly 
complex in the polaronic regime (7 = 2.5t). The struc- 
ture of <r'(u}) is different from the optical conductivity 
expected for a small polaron in the adiabatic regimeEZI. 



In particular, no dominant peak is visible at the energy 
corresponding to the depth of the lattice potential which 
traps the electron (w « 2j 2 /il). This discrepancy is 
probably the consequence of the non-adiabatic phonon 
frequency Q — t used in our calculations. In this regime, 
many-phonon optical excitations become as important 
as the purely electronic transition at u> « 2"f 2 /il. For 
weak electron-phonon coupling (7 = 0.8i) the conductiv- 
ity a' (uj) has very little weight and its structure is mostly 
determined by the discrete electronic energy levels of the 
non-interacting system.c3 

Our results agree perfectly with the known features 
of the single-electron Holstein model discussed at the 
beginning of this section. In particular, they confirm 
that the ground state is formed by an itinerant quasi- 
particle (quasi-free electron or polaron) for all electron- 
phonon couplings, but that the crossover to the polaronic 
regime is accompanied by a substantial enhancement of 
the quasi-particle effective mass. In the polaronic regime 
our results also show that the weight of incoherent pro- 
cesses is much more important than the Drude weight 
in the optical conductivity. Therefore, small perturba- 
tions which are neglected in the Holstein model, such as 
disorder or thermal phonons, are likely to suppress the 
coherent motion of the small polaron in more realistic 
models or actual materials. 



B. Two electrons with opposite spins 

The case of two electrons of opposite spins-is not as well 
understood as the single electron case.0i3'E3 For strong 
enough coupling electrons are trapped by the lattice dis- 
tortion that they generate and form an itinerant bound 
pair called a bipolaron. If both electrons are localized 
on a single site, the bipolaron is said to be small. It 
is known that a small bipolaron is formed in the strong 
coupling limit of the Holstein model with two electrons. 
What happens at weaker coupling before the onset of 
small bipolaron formation is still debated. Three differ- 
ent scenarios are possible: both electrons remain free, 
two independent polarons are formed, or a "large" bipo- 
laron is formed (a bipolaron whose spatial extension is 
larger than one site). 

Figure [| shows the most important optimal phonon 
wave functions for each possible electronic occupation of 
a site (N a = 0,1,2) for different electron-phonon cou- 
plings. (In this system, optimal phonon states are simi- 
lar if the site is occupied by one electron with spin up 
or down, so we do not distinguish between these two 
cases). As in the single-electron system, the optimal 
phonon states are just the ground state of the operator 
b'b in the weak coupling limit. The weights of these op- 
timal states are determined by the probability of a given 
site being empty (Wo = 25/36), occupied by one electron 
(W 1 = 10/36), or by both electrons (W 2 = 1/36) when 
there are two independent electrons on a six-site lattice. 
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FIG. 4. Two-electron system: most important optimal 
phonon modes as a function of the oscillator coordinate q 
and their weights W for the three possible occupations of a 
site (N a = 0, 1, 2) and different electron-phonon couplings 7. 
For a singly occupied site (N a = 1) the optimal states are 
identical for an electron with spin up or down. The weight 
Wi is the sum of the probability of the most important state 
for each of the two spin orientations. 

As 7 increases, W\ decreases and tends to zero in the 
strong coupling limit while both Wq and W2 increase and 
tend to 5/6 and 1/6, respectively. These weights corre- 
spond to the occupation probability of a particular site 
when the two electrons form an itinerant tightly bound 
pair, as in a small bipolaron. For intermediate couplings, 
the combined weight of the most probable phonon states 
in Fig. U is smaller than 1, showing the importance of 
higher optimal phonon states, but always remains larger 
than 98% of the total weight. The wave functions for 
N a = 0, 2 have large overlaps with the ground state of a 
shifted harmonic oscillator for any coupling. For 7 > 1.5t 
this shift is very close to the strong-coupling theory pre- 
diction 2N a j/£l. Therefore, these optimal states can be 
obtained with the Lang-Firsov transformation ( [l2| ) using 
an appropriate parameter g, but g reaches the theoreti- 
cal value only at strong coupling, as in the single- 
electron case. The optimal wave function for N a = 1, 
however, changes significantly as the coupling increases. 
In Fig. [| one can see that this wave function is ap- 
proximately the superposition of the wave functions for 
N a = and N a = 2 when 7 = 1.5t. This wave function 
cannot be obtained by applying the Lang-Firsov transfor- 
mation to a bare phonon state. The shape of this optimal 
state can be understood as a retardation effect. Most of 
the time a site is empty or occupied by both electrons 
and the phonon modes are in the corresponding optimal 
states. The electronic states with N a = 1 are essential 
intermediate states for allowing the coherent motion of 
the bipolaron, as the Holstein Hamiltonian ([!]) contains 
only one-electron hopping terms, but they have very low 
probability. Electrons do not spend enough time in these 
states for the phonon modes to adapt. Therefore, opti- 
mal phonon states for N a = 1 are determined by the 
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FIG. 5. Two-electron system: (a) the spectral functions 
A(p = 0, u)) (thick line) and A(p = 0,lj) (thin line) for differ- 
ent electron-phonon couplings 7. (b) Total spectral weight 
n a (j) = 0) (filled circle) and n CT (p = 0) (open circle) as a 
function of the electron-phonon coupling. 

optimal phonon states for N a = and N a = 2. 

The single-particle spectral functions A(p = 0,oj) and 
A(p = 0, oj) are similar to those of the single-electron 
system at weak-coupling. A single peak is observed at 
lo w 24 because both electrons occupy the lowest one- 
particle eigenstates with energy —24, and the ground 
state wave function is 

W « 4 =OiT 4 =o J0) « gt =0iT ct =0 JO) . (26) 

In Fig. ||(a) we show both spectral functions for different 
electron-phonon couplings. The results for bare electrons 
are qualitatively similar to those observed in the single- 
electron system. As the coupling strength increases, the 
weight of the dominate peak shifts to an increasing num- 
ber of satellites peaks until no well-defined quasi-particle 
peak can be observed. Dressing the fermion operator sim- 
plifies the spectral weight structure a bit, but important 
incoherent contributions are still observable. Moreover, 
the total spectral weight becomes small at strong cou- 
pling for both bare and dressed electrons [Fig. ||(b)]. 
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FIG. 6. Two-electron system: (a) pair spectral functions 
P(p = 0, u>) (thick line) and P(p = Q,u) (thin line) for dif- 
ferent electron-phonon couplings 7. For 7 = Lit and 1.5t, 
P(p = 0, w) has been expanded by a factor 5. (b) Total pair 
spectral weight d(p = 0) (filled circle) and d(p = 0) (open 
circle) as a function of the electron-phonon coupling. 

Therefore, the ground state can no longer be described 
by two independent electrons or quasi-particles as was 
done above for weak coupling. 

At weak electron-phonon coupling the pair spectral 
functions P(p = 0, u>) and P(p = 0, u>) have a single peak 
at u> w At on a finite system, but the weight of this peak 
vanishes as 1/N in the thermodynamic limit. Figure |^(a) 
shows both functions for stronger couplings. The evolu- 
tion of the structure of P(p = 0, ui) for increasing cou- 
pling is qualitatively similar to that of the single-particle 
spectral function. The weight of the dominant peak is 
progressively shifted to satellite peaks until no dominant 
peak can be observed. However, the total spectral weight 
d{p = 0) remains almost constant (close to l/N=l/6) for 
all couplings [Fig. ||(b)]. As this value is a finite size ef- 
fect for 7 = 0, we think that the pair spectral function 
vanishes uniformly in an infinite Holstein lattice. The 
dressed pair spectral function is much more interesting 
[see Fig. |](a)]. It shows a single dominant peak for all 



electron-phonon couplings, though small satellite peaks 
are present at intermediate coupling. Moreover, the total 
spectral weight d(p — 0) increases with 7 and tends to 
1 for strong coupling [Fig. ^(b)]. This means that the 
ground state is given by 

|^o) « At =0 |0) (27) 

in the strong-coupling regime. This state describes an 
itinerant small bipolaron with momentum p — 0. As 
the optimal phonon states for N a — 0, 2 are close to 
the states generated with the Lang-Firsov transforma- 
tion dl2| ) for strong-enough coupling, the dressed pair 
operator A p is similar to the "bipolaron operator" built 
with this transformation in Ref. 4. The structure of both 
P(p = 0, ui) and P(p — 0, ui) can be understood using 
the same arguments as for A(p — 0, ui) and A(p = 0, ui) 
in the single-electron case. Therefore, P(p = 0,u>) and 
P(p = 0, ui) contain only peaks with spacing f2 (with 
some peaks to small to be seen) starting from — Eq, where 
Eq is now the ground state energy of the two-electron sys- 
tem. For 7 = 1.5t, the position of the dominant peak of 
P(p = 0, ui) shown in Fig. ||(a) gives a bipolaron energy 
close to the strong-coupling result E = — A-f 2 /n. 

Figure |^(a) shows the evolution of the Drude weight 
and of the kinetic energy per site as a function of the 
electron-phonon coupling. In this case too, the units are 
chosen so that both quantities are equal in the absence 
of incoherent contributions to the optical conductivity. 
Both quantities decrease smoothly as the coupling in- 
creases. As in the single-electron case, the small reduc- 
tion of D at weak coupling shows the slightly renormal- 
ized effective mass of the quasi-free electrons while the 
small but finite value of D at strong coupling shows 
that the bipolaron is a heavy itinerant quasi-particlc. 
The diminishing ratio D/T means that incoherent pro- 
cesses become more important as 7 increases. The in- 
coherent part of the optical conductivity cr'(u>) is shown 
in Fig. 0(b) for both the quasi-free-electron regime and 
the bipolaronic regime. In the quasi-free-electron regime 
(7 = 0.3t), cr'(w) has a very low weight but a simple 
structure which is determined by the discrete electronic 
levels of a non-interacting six-site lattice. In the bipola- 
ronic regime (7 = 1.5£), o'{ui) is fairly complex and the 
features predicted for a bipolaron in the adiabatic regime 
are not visible.EZI For instance, there is no clear peak at 
the energy u> w A-f 2 /tt corresponding to the depth of the 
lattice potential trapping the electrons. As already sug- 
gested in the single-electron case, this is probably due to 
the non-adiabatic phonon frequency (57 = t) used in this 
work. 

Our results show that the ground state of the two- 
electron system is composed of two independent quasi- 
free electrons at least up to 7 = 0.8t but it is a small 
bipolaron at least from 7 = 1.3i. There is a smooth 
crossover from one regime to the other as the electron- 
phonon coupling increases. The nature of the ground 
state for 0.8i < 7 < 1.3i is not directly determined by 
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FIG. 7. Two-electron system: (a) the Drude weight D 
and the kinetic energy per site T as a function of the elec- 
tron-phonon coupling 7. (b) The incoherent part of the 
optical conductivity cr'(u>) in the quasi- free electron regime 
(7 = 0.3£) and in the bipolaronic regime (7 = 1.5t). 

our study. Clearly, there is no sign of a polaron in any 
of the quantities we have calculated and the small po- 
laron appears at stronger coupling (7 w 1.5£) than the 
small bipolaron. Thus, we can conclude that a pair of 
small polarons would be instable with respect to the for- 
mation of a small bipolaron. The only reasonable can- 
didate for the ground state in the intermediate regime 
is the large bipolaron, but this state is difficult to study 
in the small lattice considered here. DMRG calculations 
for long Holstein chains confirm that the ground state 
evolves smoothly from a pair of quasi-free electrons to 
a large bipolaron and then to a small bipolaron as the 
electron-phonon coupling increases .E3 As in the polaron 
case, the large effective mass of the bipolaron and the 
dominance of incoherent processes in its optical conduc- 
tivity mean that coherent motion of a small bipolaron is 
unlikely in more realistic models. 



C. Half filling 

Mean-field theory predicts that the ground state of the 
half-filled Holstein model is a doubly degenerate CDW 



state with a dimerized lattice and a gap at the Fermi 
surface for any finite electron-phppon coupling due to 
the well-known Peierls instability.El This result is exact 
in the adiabatic limit r2 — ^ and some early works 
suggested that the ground state was also a Peierls in- 
sulating state, for arbitrary coupling at finite phonon 
frequency.EZfEj More recent results suggest that quan- 
tum lattice fluctuations can destroy the Peierls insulaU 
ing state at weak coupling or large phonon frequency£3't£l 
Furthermore, it is known rigorously that the ground state 
of a finite half-filled Holstein lattice is not degenerate 
for any finite electron-phonon coupling and any finite 
phonon frequencytil Therefore, the ground state of the 
Holstein model can be a true Peierls state with CDW or- 
der and lattice dimerization only in the thermodynamic 
limit. We also note that, in the Holstein model of spinless 
electrons at half filling, it is well established that a metal- 
insulator transition occurs at finite values.of. the .aLectron- 
phonon coupling and phonon frequencyfl'LjEEro'LjO We 
have carried out a study of the ground state of the one- 
dimensional Holstein model for electrons with spin-i at 
half filling using the DMRG method. Calculations of 
static correlation functions in long chains (up to 80 sites) 
clearly show that there is a transition from a metallic 
ground state to a Peierls insulating state with long range 
CDW order in the thermodynamic limit. For f2 = t this 
transition occurs between 7 = 0.8t and 7 = 0.9t. In this 
paper we concentrate on the dynamical properties of the 
six-site system and our DMRG results will be reported 
elsewhere. til .— ■ 

We had previously noticecO that static correlation 
functions reveal a crossover from a uniform ground state 
to a Peierls CDW ground state in a half-filled six-site 
Holstein lattice around 7 = f for £1 = t. In fact, despite 
the absence of a true broken symmetry ground state in a 
finite Holstein lattice, signs of this crossover and the ex- 
istence of an "insulating" Peierls CDW phase are clearly 
seen in static and dynamical properties. This is due to 
a quasi-degeneracy of the ground state at strong enough 
coupling. In Fig. ^(a) we show the energy difference be- 
tween the ground state |^o) and the first excited state 
\tpi). Above 7 = Lit the difference is very small and the 
two states are almost degenerate. The eigenstates | Vo) 
and \ipi) have momentum and n, respectively. We know 
that they have a constant density (V'o.i l n i|"0o,i) = 1 an< ^ 
a uniform lattice structure (Vfyil&l + &i|^o,i) — 27/fi, be- 
cause of the translation symmetry of the Holstein Hamil- 
tonian (|l|) with periodic boundary conditions. If they 
were exactly degenerate, we could build two broken sym- 
metry eigenstates \ip±) — \tpo) ± \ifji) with a charge mod- 
ulation (ip±\rii\ip±) = 1 ± (— iy6 and a dimerized lat- 
tice (tfj ± \b\ + bi\^ ± ) = 27/fi (1 ± (-1)^). These two 
states would correspond to the two possible phases, of 
the degenerate Peierls CDW state obtained in mean-field 
approximation or in the adiabatic limit .113 As the two- 
lowest eigenstates \ipo) and are quasi-degenerate for 
7 > Lit, the six-site system properties are almost indis- 
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FIG. 8. Half-filled system: (a) gap between the ground 
state energy E$ and the energy E\ of the first excited state 
and (b) density of doubly occupied sites d as a function of the 
electron-phonon coupling 7. 

tinguishable from those of a true degenerate Peierls CDW 
state. 

The optimal phonon wave functions for the half-filled 
Holstein model have already been discussed in our pre- 
vious work.E3 As expected, they are similar to the wave 
function obtained in the two-electron system and only 
the relative weights of the optimal states are quite dif- 
ferent due to the different number of electrons in the 
system. The total weight of optimal phonon states as- 
sociated with zero or two electrons on a site tends to 1, 
while the weight of those associated with one electron on 
a site becomes very small for strong coupling. This aspect 
of the crossover from independent electrons to the Peierls 
CDW state is also illustrated by Fig. ||(b). The density of 
doubly occupied electronic sites increases rather sharply 
from the independent electron result 1 /4 to the maximal 
value 1/2 possible in a half-filled system. Thus, above the 
crossover regime the electrons form tightly bound pairs. 
These electronic pairs are heavily dressecLby phonons and 
thus, can be seen as small bipolarons.Ej In the strong- 
coupling anti-adiabatic limit, it is known exactly that at 
half-filling these small bipolarons form an ordered phase, 
which isJjully equivalent to the Peierls CDW state in this 
regime. Ell 

The spectral weight functions of the half-filled system 
are harder to interpret than those of the single- and two 
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FIG. 9. Half-filled system: (a) the spectral functions 
A(p — 7r/3) (thick line) and A(p = 7r/3) (thin line) for 
different electron-phonon couplings 7. (b) The momentum 
density distribution n a (p) for all values of the momentum p 
(results for p and — p are identical) as a function of the elec- 
tron-phonon coupling. 

-electron system. Nevertheless, some of their proper- 
ties can be understood. In the non-interacting limit, the 
spectral functions A(p, ui) and A(p, to) show a single peak 
with total weight 1 for momentum \p\ < 7r/3 correspond- 
ing to occupied single-electron states and are uniformly 
zero for momentum \p\ > 7r/3 corresponding to unoccu- 
pied single-electron states. As the electron-phonon cou- 
pling increases, both A(p, ui) and A(p, ui) become fairly 
complex for all values of p. As an example, Figure |^(a) 
shows both functions for p = 7t/3. In the strong-coupling 
regime (7 = 1.5t), one can see that the highest peak of 
the bare spectral function A(p, u>) is located at an en- 
ergy ui w 47 2 /f2 corresponding to the energy required 
to remove one electron of the half-filled system without 
perturbing the dimerized lattice structure of the Peierls 
CDW ground state. On the other hand, the spectral 
weight for dressed electrons A(j>, ui) shows a relatively 
sharp dominant peak, but with only a fraction of the 
total weight, for all couplings. The position of this peak 
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FIG. 10. Half-filled system: (a) the pair spectral functions 
P(p = 0,w) (thick line) and P(p = 0,w) (thin line) for dif- 
ferent electron-phonon couplings 7. (b) The total spectral 
weight d(p) for all values of the momentum p as a function of 
the coupling 7. 

corresponds to the difference between the ground state 
energy of the half-filled system and that of a system 
with one electron removed from a half-filled band. For 
7 = 1.5£ we find that this peak is around u> = 3j 2 /Q in 
agreement with strong-coupling theory predictions. The 
total spectral weight n a (p) is shown in Fig. ||(b) as a 
function of 7. Results for n a (p) are similar. One can see 
that n a (p) evolves smoothly from the momentum density 
distribution of free electrons, n a (p) = 1 for \p\ < 7r/3 and 
n <j(p) = for \p\ > 7r/3, to that of completely localized 
electrons, n a {p) — 0.5 for all p, a. 

At weak coupling, the properties of the pair spectral 
functions P(p, w) and P(p, to) depend strongly on the 



momentum p but can be easily understood from weak- 
coupling perturbation theory. In the Pcicrls CDW regime 
both functions acquire a simple structure, which is largely 
similar for all values of the momentum p. For instance, 
Figure |T^(a) shows both functions for p — from weak 
to strong coupling. For very weak coupling the bare 
pair spectral function P(p = 0, to) has only two peaks 
at lo — 2t and At, which correspond to zero momentum 
pairs made of two electrons with momentum p = ±7r/3 
and p = 0, respectively. As the coupling increases toward 
the crossover regime, the spectral weight is shifted to an 
increasing number of satellite peaks. In the Peierls CDW 
regime we observe a relatively broad cluster of peaks. 
This cluster seems to be centered around an energy which 
is slightly lower than the energy 8j 2 /fl required to re- 
move two bare electrons from the half-filled ground state 
without disturbing the dimerized lattice structure. For 
all couplings P(p — 0,lu) has a well-defined dominant 
peak, although incoherent contributions are significant in 
the quasi-free-electron regime. The position of this peak 
gives the energy difference between the half-filled band 
ground state and the ground state of a system with two 
electrons removed from a half-filled band. For 7 = 1.5t 
one can see that this energy difference is about 4j 2 /£l in 
agreement with strong-coupling theory. As noted above 
similar results are obtained for all other momentum in 
the Peierls CDW phase. 

The well-defined peaks observed in P(p, to) confirm 
that the ground state is made of small bipolarons in this 
regime. In such a case, the total weight d(p) can be seen 
as the bipolaron momentum distribution. The distribu- 
tion d(p) « 0.5 for all momentum p that we observe at 
strong-coupling [see Fig. [l0|(b)] indicate that bipolarons 
are completely localized in the Peierls CDW regime. The 
study of static correlation functions shows that these lo- 
calized small bipolarons form an ordered phase even for 
finite electron-phonon coupling and finite phonon fre- 
quencyJ13 Therefore, we think that the ground state of 
the half-filled six-site Holstein lattice can be seen either as 
a Peierls CDW state with lattice dimerization or as an or- 
dered phase of localized small bipolarons. The first point 
of view corresponds to the adiabatic limit (f2/t — » 0) 
result and the second is more appropriate in the strong- 
coupling or anti-adiabatic limit (i/fi — > 0). For the inter- 
mediate case n = t discussed here, both pictures appear 
completely equivalent . 

Figure [ll](a) shows the evolution of the Drude weight 
and of the kinetic energy per site as the electron-phonon 
coupling increases. Again the units are chosen so that 
both quantities are equal in the absence of incoherent 
processes in the optical conductivity. The kinetic en- 
ergy decreases rather smoothly as in the single- and two- 
electron cases. On the other hand, we can clearly see a 
sharp decrease of the Drude weight around 7 = t. Above 
7 = 1.5£, D vanishes within numerical errors. (Of course, 
D is never really zero on a finite cluster because of quan- 
tum tunneling.) This negligible value of D is consistent 
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FIG. 11. Half-filled system: (a) the Drude weight D and 
the kinetic energy per site T as a function of the elec- 
tron-phonon coupling 7. (b) The incoherent part of the 
optical conductivity cr'(u)) in the quasi- free-electron regime 
(7 = 0.3t) and in the Peierls CDW regime (7 = 1.5t). 

with the ground state being made up of localized small 
bipolarons. Although a six-site cluster can not be metal- 
lic or insulating, the behavior of D shown here illustrates 
perfectly the metal-iiisiilator transition observed in the 
thermodynamic limit Ej The incoherent part of the opti- 
cal conductivity is shown in Fig. [ll](b). For weak cou- 
pling (7 = 0.3i), c'(u>) has very little weight and contains 
a single dominant peak that can be explained by the dis- 
crete electronic energy levels of the non-interacting sys- 
tem. In the Peierls CDW regime (7 = 1.5i), the structure 
of cr'(u>) is more complex but clearly shows a significant 
peak around an energy oj « A^ 2 /Q corresponding to the 
Peierls gap in the strong-coupling limit. In the adia- 
batic limit the conductivity is zero below the Peierls gap. 
The significant tail observed below the peak in Fig. O(b) 
is due to phonons. It is an excellent illustration of the 
subgap optical absorption predicted in-JPeierls systems 
because of quantum lattice fluctuations!^ 

Our results clearly show that there are two distinct 
regimes in the half-filled six-site Holstein lattice. The 
crossover from one regime to the other occurs rather 
sharply at a critical electron-phonon coupling, which is 



7 w t for = t. Below this critical coupling, the 
static and dynamical properties of the system are those of 
quasi-free electrons on a finite cluster. Above the critical 
coupling, the ground state is made of ordered localized 
small bipolarons, which can also be seen as the Peierls 
state with CDW order and lattice dimerization predicted 
by mean-filed theory. These two regimes are the finite 
system precursors of the metallic and insulating Peierls 
CDW phase of the infinite system, and the crossover be- 
tween both regimes signals the quantum phase transition 
observed at finite electron-phonon coupling in tha-ther- 
modynamic limit of the half-filled Holstein modeLE^ 



IV. CONCLUSION 

We have studied the dynamical properties of the six- 
site Holstein model for three different electron concentra- 
tions using an exact diagonalization technique. A density 
matrix approach allows us to generate an optimal phonon 
basis and to truncate the phonon Hilbert space with- 
out significant loss of accuracy. With this very efficient 
method we are able to observe the evolution of the system 
properties as one goes from the weak electron-phonon 
coupling regime to the strong-coupling regime. For all 
three electron concentrations studied a smooth crossover 
is observed from quasi-free electrons at weak coupling to 
a strongly correlated state in the strong-coupling regime. 
This strongly correlated state is a heavy itinerant small 
polaron in the single-electron heavy itinerant 

small bipolaron in the two-electron case, and a set of or- 
dered localized small bipolarons similar to a Peierls CDW 
state in the half-filled band case. 

The study of the optimal phonon states reveal that 
they often are the eigenstates of a quantum oscillator 
with a shifted equilibrium position. These states can 
be obtained by applying the Lang-Firsov transformation 
to the bare phonon states. The amplitude of the shift 
is generally smaller than the exact results for t = but 
tends to this value for strong coupling. However, in some 
cases we have observed some retardation effects due to 
the electron motion and the slow response of the phonon 
modes. These effects are small but they are essential 
for an accurate description of the electron motion along 
the lattice, either as independent particles or as part of 
a composite quasi-particle (polaron or bipolaron). Al- 
though we have presented only results for £1 = t, we have 
checked that these retardation effects become more im- 
portant for smaller values of £l/t (adiabatic limit) but 
vanish for larger values of Vt/t (anti-adiabatic limit). 

We have obtained a wealth of information from the 
single-particle and pair spectral weight functions. By 
dressing the electron operators with the optimal phonon 
states, we are often able to simplify the structure of 
these spectral functions and to obtain well-defined quasi- 
particle peaks. For instance, we can identify the for- 
mation of a small polaron and of a small bipolaron by 
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the appearance of a single dominant peak in the single 
electron spectral function A(u>) and in the pair spec- 
tral function P(u), respectively. We have also studied 
the optical conductivity for these systems. In all cases 
the Drude weight decreases substantially as the electron- 
phonon coupling increases from the non-interacting limit 
to the strong-coupling limit. At half filling, the Drude 
weight decreases abruptly around 7 = t and is negligibly 
small for larger 7. This suppression of coherent transport 
is linked to the appearance of a quasi-degenerate Peierls 
CDW ground state. These results support measurements 
of static correlation functions in larger systems which 
shows the existence of a metal-insulator transition in the 
thermodynamic limit of the half-filled one-dimensional 
Holstein model. 

One obvious limitation of this study is the complete ne- 
glect of electron-electron interaction when there is more 
than one electron on the lattice. This interaction is likely 
to strongly affect the properties of the system, especially 
in the regime where bipolarons are formed in the ab- 
sence of electronic repulsion. The approach used in this 
work can be applied without difficulty to models with 
an electron-electron interaction and we have started to 
investigate the Holstein-Hubbard model, which includes 
an on-site electron-electron repulsion. 
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